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Abstract 

The scale dependence of the ratios of parton distributions in a proton of a nucleus A and 
in the free proton, Rf{x,Q'^) = fi/A{x, Q"^) / fi{x, Q^) , is studied within the framework of 
the lowest order leading-twist DGLAP evolution. By evolving the initial nuclear distribu- 
tions obtained with the GRV-LO and CTEQ4L sets at a scale Q^, we show that the ratios 
Rf{x, Q^) are only moderately sensitive to the choice of a specific modern set of free parton 
distributions. We propose that to a good first approximation, this parton distribution set- 
dependence of the nuclear ratios Rf{x, Q"^) can be neglected in practical applications. With 
this result, we offer a numerical parametrization of Rf{x, Q^) for all parton fiavours i in any 
A > 2, and at any 10~^ ^ x < 1 and any > 2.25 GeV^ for computing cross sections of 
hard processes in nuclear collisions. 
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1 Introduction 



The measurements of the nuclear structure function Fc^i^x^Q"^) in deeply inelastic lepton- 
nucleus scattering (DIS) indicate clearly that parton distributions of the bound protons 

are different from those of the free protons, fi/A{x,Q^) ^ fi/p^x^Q"^). The nuclear effects 
are often categorized according to those observed in the ratio of the structure functions 
of nuclei relative to deuterium, Rp^ = F^/Fj^: shadowing {Rp^ < 1) at Bjorken-x < 0.1, 
anti-shadowing (i?^ > 1) at 0.1<x<0.3, EMC effect (i?^ < 1) at 0.3<x<0.7, and Fermi 
motion {Rp^ ^ 1) towards x — > 1 and beyond. The recent high-precision measurements by 
the New Muon Collaboration (NMC) of the structure function F2 of tin vs. that of carbon, 
Fi"^ I F2 [H, have also revealed a Q^-dependence at small values of x. 

Theoretically, the origin of the nuclear effects is still under debate but it is believed 
that different mechanisms are responsible for the effects in the different regions of x. For a 
compact introduction and references to the various theoretical models we refer the reader to 
Arneodo's review, Ref. 

In this paper, as a sequel to Ref. |jl0[, we are going to focus on studying the perturbative 
QCD-evolution of the ratios Rf = fi/A{x,Q^)/ fi/p{x,Q^) within the framework of lowest 
order (LO) leading twist (LT) DGLAP evolution [0]. To our knowledge the perturbative 
QCD-evolution of nuclear parton densities has been studied at least in Refs. |[T2|-||22||. 

Our approach is in practice very similar to the case of the free proton: once the absolute 
nuclear parton distributions are given at an initial scale = Ql ^ ^qcd and at some 
^^rnin ^ X < 1, the further evolution in is predicted (at this range of x) by the DGLAP 
equations. In other words, as in |jlO[, we assume that at Q^^Ql the scale evolution of the 



nuclear parton densities is purely perturbative and that the initial conditions at Qq contain 
nonperturbative input. We determine the initial nuclear parton distributions iteratively 
through the QCD evolution by using experimental data and conservation of momentum and 
baryon number as constraints. It should be noted that we do not try to explain the origin 
of the nuclear effects but we study their behaviour once they are there at an initial scale. 
Furthermore, since we neglect all the higher twist contributions, like the perturbative GLR- 
MQ terms |j23| at small values of x, no nuclear effects are generated through the evolution 
but they are all hidden in the initial distributions at Qq. 

In Ref. [|ig, the DIS data on F^/F^ and F^/F^ 0-0 were used together with the 



measurements of the Drell-Yan (DY) cross sections in pA vs. pD collisions and with the 
conservation of baryon number and total momentum to determine the initial nuclear parton 
distributions at a scale Ql = 2.25 GeV^. Especially, it was shown in [|lOl that the LT LO 
DGLAP-evolution can account very well for the dependence observed in the experimental 
ratio Fi'^/F^ g. 

The measurements of the ratios Rp^ in DIS, and also the ratios of nuclear Drell-Yan 
(DY) cross sections relative to deuterium, provide us with information on the nuclear parton 
distributions only relative to the nucleons in deuterium. In an approximation where the small 
nuclear effects in deuterium are neglected, as done in |Ty] and as we will do here as well, we get 
information relative to the free nucleons. Then, quite obviously the extraction of the absolute 
nuclear parton distributions depend on our choice for the set of parton distributions in the 
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free proton, and the QCD-evolution itself will depend on this choice, tooQ. In Ref. 
it was anticipated that the ratios Rfi^x^Q"^) should not depend as strongly on the choice 
for the set as the absolute distributions do. In this paper, we will verify this statement 
quantitatively by comparing the scale-evolved nuclear effects initially obtained by using the 
GRV-LO [H and CTEQ4L ^ distributions. 

After demonstrating that the set-dependence in the nuclear ratios Rf{x, Q"^) for individ- 
ual flavours indeed is a negligible effect as compared to the current overall uncertainties, it 
becomes more meaningful to present a numerical parametrization of Rf{x, Q"^) for practical 
applications. The parametrization can then be used with any modern set of LO parton 
distributions to obtain the absolute nuclear parton distributions needed for computing hard 
scattering cross sections in nuclear collisions. The parametrization we offer can be called for 
any parton flavour any A > 2, 10~^ < x < 1 and > 2.25 GeV^. The parametrization, 
which we give in form of a short Fortran code, is now available from usQ. 



2 The framework and the assumptions 



A detailed formulation of our approach can be found in Ref. ||T0[ but let us recall the un- 
derlying basic assumptions also here. The starting point is the measured ratio of structure 
functions of a nucleus A and deuterium D, Rp^ = F.f/F^. In the one-photon approximation 
the ratio of the cross sections measured in I A and in ID gives directly the ratio Rp^^ P™" 
vided that the ratio of the cross sections of longitudinally and transversally polarized virtual 
photons with the target, R = ai/cTT, does not depend on the target. So far, experiments 
have not shown any significant target dependence ^, |[. Our approach, being currently 
only in the lowest order in the cross sections, will not generate such a dependence, either. 

First, the structure function ratio Rp^ = F.^/F^ is expressed in terms of nuclear par- 
ton distributions by using the LO relation in the QCD improved parton model, F2 = 
J2q ^1[xfg{x, Q'^)+xfg{x, Q^)]. By the nuclear distributions fi/A of a parton flavour i we mean 
the average distributions of a parton flavour i in a proton of a nucleus A: Q^) = 

fi/Ai^jQ"^). In principle, the nuclear parton distributions are non-zero up to x — >■ A but 
the small tails are completely negligible from the point of view of our study. We therefore 
approximate fi/A{x > l,Q^) = 0. 

For isoscalar nuclei dn/A = ^p/A = ua and Un/A = dp/A = dA obviously hold but we have 
to assume that these are good approximations for non-isoscalar nuclei as well^. We also 
neglect the small nuclear effects in deuterium |^. With these approximations Rp^ reflects 
directly the deviations of F^ from the free nucleons. 

The nuclear effects in the parton distributions are defined through the ratios for each 
flavour as, 

fi/A{x,Q^) 



Rtix.Q' 



^We will refer to the dependence on the choice for the set of parton distributions in the free proton as 
"set-dependence" 



^via email or in the WWW, http: / /fpaxpl. usc.es/phcnom/ 

^For non-isoscalar nuclei the corresponding exact relation is d{u)n/Az — '^{.'^pIA^ 
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where is the distribution of the flavour i in the free proton as given by the chosen set of 
parton distributions. The modern sets include the rapid increase of the gluon and sea quark 
distributions at small values of x, in accordance with the rise discovered in the structure 
function at HERA p8[. In the DGLAP evolution we treat the massive quarks as massless 



and generate them only radiatively above fixed threshold scales. The two modern LO sets of 
parton distributions in which the heavy quarks are treated in this manner are the GRV-LO 
| 25[| and CTEQ4L ||26|, which is why we use only these two sets in this work. 

To eliminate the complication of determining nuclear effects for the heavy quark distri- 
butions at the initial scale Ql, we choose Ql at or below the charm-threshold Q^, where 
xc = xcA = 0. In the set GRV-LO the threshold is Ql = 2.25 GeV^ and in the CTEQ4L 

= 2.56 GeV^. To be consistent, we use these in the corresponding evolution in the 
nuclear well. The HERA results have also shown that the leading twist evolution 

describes the structure function well at Q^^l GeV^ and 10~^ p9[. In the nucleus, 



however, the higher twist effects can be expected to be stronger p3| , [12| , |T^. Since we will 



not include these here, we should not start the perturbative evolution of nuclear parton 
distributions below ^1 GeV^. Therefore, for not too small values of x, it is fairly safe to 
choose Ql = Ql = 2.25 GeV^ of the set GRV-LO. The results in Ref. ^ show that this is 



a valid initial scale at least for x ^ 0.01, where the higher twist contributions are negligible. 

As further approximations, but for the initial conditions at = Qq only, it is assumed 
that the differencies between the nuclear effects for sea quarks and antiquarks, as well as for 
different sea quark flavours are negligible: R^{x,Ql) = R^^^^{x,Ql) = Rg{x,Ql). Similarly, 
the same nuclear effects are assumed for the valence quarks at the initial scale, R^^{x, Ql) = 
R^^{x, Ql) = Ry{x, Ql). The fact that with these approximations we get a stable evolution 
- i.e. that the nuclear effects do not rapidly evolve away from what is assumed at Ql - shows 
that these are indeed reasonable first approximations. 

Then, we arrive at a simple formula for an isoscalar nucleus, 

Ri{x,Ql) = A\fix,Ql)R^{x,Ql) + [Ai',{x,Ql) + Mx,Ql)]Rj{x,Ql), (2) 

where the coefficients are known exactly |[T0|| : 

A\f{x,Q^) = 5[uv{x,Q^) + dv{x,Q^)]/NF,{x,Q^) (3) 

Aii{x,Q') = 10[u{x,Q') + d{x,Q')]/NF,{x,Q') (4) 

A{x,Q') = As{x,Q')/Np,{x,Q') (5) 

Np.ix^Q^) = 5[My(a;,Q') + c?y(a;,Q')] + 10[m(x,Q2) + J(x,g2)] +4s(x,g2). (6) 

They obviously depend on the specific set of parton distributions chosen for the free proton. 
The above Eq. (|]) correlates the nuclear effects of the valence quarks with those of the sea 
quarks through the ratio Rp^. As described in detail in ||T0[, we fix Rp^ iteratively through 



the evolution by using the DIS data g-f^, decompose i?^ i^^^ and R^ at = q2 
according to Eq. by simultaneously constraining Ry by conservation of baryon number 
and Rg by the DY data [^. At large values of x the sea quark distributions are very 
small as compared to the valence distributions, and the contribution from the sea to Rp^ is 
negligible. This is why we have no constraints for Rg at x^ 0.3. We simply assume that the 
sea quarks show a similar EMC effect as the valence quarks do, so at large values of x we 
have R^ = R^ = i?^. 



3 



Also the region below x^O.Ol is not well constrained: there exist experimental DIS 
data, even at much lower values of x 0, ||, but due to the correlation of x and (Q^) in the 
experiments, the values of (Q^) are in the nonperturbative region from our point of view. 
Also, as pointed out in 0, since dR^^{x, Q^)/d\ogQ'^ ~ RqC^x, Q^)-Rj^^{x, Q^), the slope 
of Rp^ in - even the initial sign of the slope - depends on what is assumed for the initial 
nuclear gluon profile Rq{x, Qq) at small values of x. We assume a saturation of shadowing in 
Rp^, as observed in the nonperturbative region P, P], and again, to reach a stable evolution 



Tod, take Rq = Rp^ at small values of x. In relation to the data on Rp^ measured in the 



region of very small x |^, we basically expect that the nonperturbative evolution does not 
change the sign of the Q^-slope of Rp^. This interesting question is related to the origin of 
nuclear shadowing. Within our approach, however, we cannot address this problem further. 

Apart from the overall momentum conservation, we do not, unfortunately, have any 
direct constraints for the nuclear gluon distributions. In principle, as pointed out in 
the dependence of Rp^, if measured very accurately at small enough values of x, would 
serve as such a constraint. As a step towards this direction, the NMC has recently measured 
the Q2 dependence of Fi^'/F^ g and Pirner and Gousset |^ have used this data to 
extract the approximate nuclear gluon densities. Since we now assume a similar shadowing 
in Rq as in Rp^, we are bound to have some antishadowing due to the overall momentum 
conservation ||Ty]. In our analysis, we use the results of Ref. |^, to constrain only the point 
where Rq{x,QI) ~ 1. By requiring again a stable evolution (unlike in [jl4|) we expect that 
the gluons have an EMC-effect as well. The height of the antishadowing peak in Rq{x, Qq) 
is then fixed by the momentum conservation. As shown in |T^, with the full scale evolution 
this results in quite a good overall agreement with the analysis in Ref. PU|. 



3 The set-dependence of the nuclear effects 

In this section, we will make a quantitative comparison of the set-dependence of the nuclear 
effects for all parton flavours. More exactly, we will compare the results of [|l^ obtained 



with the parton distribution set GRV-LO pH| to those obtained with the set CTEQ4L 



First, the set-dependence will obviously be the largest for largest nuclei, so it is sufficient 
to focus on the set-dependence for an isoscalar A = 208. To have a consistent but the 
most straightforward comparison, we fix the parametrization of Rp^{x, Ql) to be exactly the 
same as given in ||10|. By doing this, we anticipate that the set-dependence of the scale 



evolution of Rp^ will be negligible, so that a new iteration with a new initial parametrization 
of Rp^{x,Ql) is not needed. Similarly, we keep the parameters for Ry{x,QI) at x<0.1 the 
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same as given m 

To apply the initial parametrization of Rp^{x,Ql) to the CTEQ4L distributions in a 
consistent manner, and to remove a source of uncertainty, we first evolve the CTEQ4L 
distributions downwards from = Qq ^"'"^'^^^ = 2.56 GeV^ to our starting scale = Ql = 
2.25 GeV^, and only then apply the procedure of extracting the initial nuclear effects in the 
valence quarks, Ry{x,Ql), in the sea quarks, Rg{x,Ql), and in the glue, Rq{x,QI). These 
initial ratios can be found in Fig. |l|, together with the previous results |]lOl with the set 
GRV-LO. The resulting set- dependence in Rg{x,Ql) and Ry{x,Ql) is small and certainly 
well within the expected overall uncertainty in fixing the initial conditions for these ratios. 
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Notice also that the set- dependence of the initial gluon ratio Rq{x,Qq) is conveniently 
small, as seen in Fig. |l|, even though the absolute gluon distributions of GRV-LO and 
CTEQ4L differ from each other, even by a factor ~2. Remember that the shadowing of 
Rq{x ^ Q o) at very small values of x was fixed to coincide with Rp^{x, Qq) but the antishad- 
owing peak in Rq{x, Qq) is determined from the overall momentum conservation. 

Next, we perform the DGLAP evolution |Tl|] for the absolute nuclear parton distributions 
now obtained by using the set CTEQ4L, and form the nuclear ratios of individual ratios at a 
large scale = 10000 GeV^. The comparison with the ratios obtained by the set GRV-LO 
is shown in Fig. |l|. It is quite interesting to notice that even if the absolute distributions may 
differ from each other considerably, in the scale-evolved nuclear ratios the differencies tend 
to cancel out, making the set-dependence a small effect. The set-dependence is the largest 
for the gluon ratios at small values of x, and this is transmitted directly to the sea quarks 
and to the heavy quarks. The strange-quarks also show a set-dependence of at most 5 %. 
Another interesting point is that the difference between R^ and Rf stays quite small even 
though M 7^ (i for the set CTEQ4L. 

Bearing in mind the uncertainties in fixing the initial ratios, especially the assumptions 
regarding shadowing at x<0.01, and the gluon profile in general, the conclusion is that 
the set-dependence is a negligible effect as compared to the current overall uncertainty in 
determining the initial conditions. 

In Fig. 0, the set-dependence of the nuclear ratios is shown as function of nuclear mass 
A at different fixed values of x and one fixed value of Q^. This figure demonstrates the fact 
that for smaller nuclei the set- dependence is even a smaller effect. 



4 The parametrization for practical applications 

Originally [|14|, [l^, we have stored the x- and Q^-dependent nuclear ratios, Rf{x, Q^) for all 
relevant parton flavours in a big table for each nucleus separately. In practical applications of 
computing cross sections of hard scatterings in nuclear collisions, these big tables have to be 
interpolated to obtain the absolute nuclear distribution^ Q^) = Rf{x,Q'^)fi{x,Q'^). 

Then, in principle, one should store ~ 200 tables to be able to get nuclear parton distribu- 
tions in an arbitrary nucleus A > 2. To avoid this complication, and to make the nuclear 
effects more easily available for any user, we have prepared a numerical Fortran code which 
parametrizes the x-, Q^- and especially the A-dependence of Rf{x, Q^) for all parton flavours 
i of an isoscalar nucleus A. 

We offer the numerical parametrization of the ratios Rf{x, Q"^) in two formats: The flrst 
one, called "EKS98" , computes the nuclear ratios for all parton flavours within one function 
call, and the second one, called "EKS98r" , computes the nuclear ratio of the parton flavour 
specifled by the user. Both versions are initialized by reading in two small tables. The codes 
and the tables are contained in a package now available from any of the authors via email. 
For an easier access, the package is also placed in the WWW0. 

Here we do not wish to go into the numerical details of making such a parametrization. 
Instead, we simply demonstrate the quality of our parametrization in Figs. (|l|) and @. As 

^Also these tables with an interpolation routine is available from us via email, 
^http:/ /fpaxpl. usc.es/phenom/ 
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can be seen in the figures, the parametrization is based on the results |]T0[ obtained with 



the set GRV-LO. Therefore the parametrization is closer to those results and the region of 
validity is 10~^ < a; < 1 and > 2.25 GeV^. For A < 2, the parametrization simply 
returns unity, in accordance with the fact that nuclear effects in deuterium were neglected. 

Now, with the results of the previous section on the approximate set-independence of 
the ratios Rf, the absolute parton distributions in protons of an arbitrary nucleus A with Z 
protons and A — Z neutrons can be obtained simply by multiplying the parton distributions 
yPDset q£ ^^y- jYiodern lowest order set by our parametrization /^^'^^^^^^ 

/^r Q') - fF/^/fi^. Q') - ^(x, Q')fr-\x, Q% (7) 

and the corresponding distributions in a bound neutron can be obtained through the assumed 
isospin symmetry (exact for isoscalars, though). Thus, our parametrization should be easy 
to use together with parton distributions from e.g. the CERN PDFLIB or other similar 
packages. 



5 The set-dependence in the physical quantities 

As the last task, we check explicitly the set-dependence in the physical quantities studied in 



Ref. |]T0|, and simultaneously we show the range of uncertainty when using our parametriza- 
tion with different sets of parton distributions in the free proton instead of the nuclear ratios 
obtained specifically for the chosen set. For brevity, let us call the latter ones as GRV-ratios 
and CTEQ-ratios. 

First, in Fig. |^, we calculate the dependence of the ratio / = Rp^/Rp^ for an 
isoscalar A = 118, corresponding to an isoscalar tin, and compare our calculation with the 
NMC data There are four curves in the figure: two curves which have been obtained by 
using the nuclear parton distributions obtained with the GRV-ratios and CTEQ-ratios, and 
two curves obtained with our numerical parametrization. The four curves practically lie on 
top of each other, so the uncertainty due to the choice of the set, and the uncertainty in 
our parametrization, are negligible as compared to the experimental (statistical) error bars. 
We have checked that we arrive in the same conclusion also for other ratios Rp^{x,Q'^)^ as 
expected due to the small differencies of the ratios in Fig. |l[ 

Another quantity we study here is the ratio of the lowest order Drell-Yan cross section in 



pA and pD and its comparison to the E772-data p4|. As presented in |T0[, the cross section 



ratio for a nucleus A with Z protons and A — Z neutrons can be written as 

_ jda^/dx2dQ^ 
np,Y[X2,^ ) = Ty~pD 



= {4[ni(n^ + d^) + u,{u^ + d^)] + [d,{d^ + u^) + d,{d^ + u^)] + As,s^ + ...}/Nr>Y 
17 

+ (— - l){A[u,{u^ - d^) + ur{4 - 4)] + [d,{d^ - u^) + Ji(4 - u^)]}/NnY (8) 
where the denominator is 

Ndy = 4[ui(m2 + d2) + ui{'^2 + di)] + [di{d2 + U2) + di{d2 + ^2)] + 4siS2 + ... (9) 
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and where we have used the notation ql = Q"^) for i = 1, 2 and q = u,d, s, .... The 

scale in the parton distributions is the invariant mass of the lepton pair. The target 
(projectile) momentum fraction is X2 {xi), and Xi = Q^/{sx2). 

In Fig. i we plot the ratio R^y for ^eC, toCa, ^^Fe and ^^jW as a function of the 
target momentum fraction X2- As in Fig. ^, we show four different calculation in each 
panel, corresponding to X2 and (Q^) of the data |Q: two of them are obtained by using 
the GRV-ratios and the CTEQ-ratios with the corresponding sets for the free proton, and 
the remaining two by using our numerical parametrization together with the GRV-LO and 
CTEQ4L parton distributions for the free proton. Carbon and calcium are isoscalar nuclei, 
A = 2Z, so for them the latter part in Eq. (|^) drops out. The remaining small difference 
in the calculated results (the big circles and the diamonds) is then the uncertainty in our 
parametrization due to the set-dependence. With iron and especially with tungsten the 
situation becomes more interesting because the non-isoscalar effects start to play a visible 
role. In the panel for tungsten, we clearly see that in the calculated quantities, there are 
separate "error bands" for the results obtained with GRV-LO and CTEQ4L distributions. 
The "width" of these bands is again the uncertainty related to using our parametrization 
instead of the "exact" GRV- or CTEQ-ratios. As compared to the differencies in the free 
proton distributions and to the experimental error bars for -R^y, the uncertainty due to 
using our parametrization instead of the "exact" nuclear ratios is small. 



6 Discussion and conclusions 



ff^^^{x, Q^), obtained with different sets may differ considerably in absolute magnitude but 
in the scale dependent the nuclear ratios, Rf{x,Q'^) = /^2''''*(x, Q^)//™''''*(s, Q^) the dif- 



By using the approach of Ref . [T^ , we have studied the leading twist lowest order DGLAP 
evolution of nuclear parton densities. The nuclear effects to the initial parton densities for 
free proton at = Ql have been determined by using the data from the deeply inelastic 
lA scatterings 0-0 and from the Drell-Yan measurements in pA collisions 1 24], and conser- 
vation of baryon number and momentum as constraints. The nuclear parton distributions, 

fPDs 

h/A 

li/A 

ferencies between the sets are much smaller. By using the GRV-LO and CTEQ4L 
distributions, we have verified that the set-dependence in Rf{x,Q'^) for individual parton 
flavours is negligible as compared to the current theoretical overall accuracy in fixing the 
initial nuclear effects at Q"^ = Qq, and also compared to the precision of the experimental 
data. 

With the result that the dependence of the nuclear ratios on the choice for the parton 
distributions for the free proton is negligible, we have prepared a numerical parametrization 
of the ratios Rf{x,Q'^) to be used for any practical applications of computing cross sec- 
tions of hard scatterings in nuclear collisions, where absolute nuclear parton distributions, 
fF/A^^^i^^Q'^) — Rf{^^Q^)fF^^^^i^yQ'^) needed. We distribute our parametrization in a 
package containing two alternative Fortran routines. The package is now available from us 
via email or from the WWW0. 

In Ref. it was shown that the pure leading twist, lowest order DGLAP evolution can 



^http://fpaxpl. usc.es/phenom/ 



7 



account very well for the observed scale dependence of the structure function ratio F2^/F2- 
Here we have demonstrated by using the GRV-LO and CTEQ4L distributions that using any 
modern LO set of parton distributions for the free proton will result in the same conclusion. 

We have also shown that quite different isospin effects can be expected in the ratios 
of Drell-Yan cross sections in pA vs. pD for large non-isoscalar nuclei like tungsten when 
using sets with different relative magnitudes for u and d. We have also demonstrated that 
it is again sufficient to use set-independent nuclear ratios, i.e. our parametrization, in the 
computation of the DY ratios, the differencies result from the differencies between the parton 
distributions of the free proton. 

We conclude by recalling the questions in our approach that are still open regarding the 
determination of the initial distributions at = Qq- 

At large values of x the existence of the EMC-effect in the sea and for the gluon dis- 
tributions is an assumption. To directly verify this experimentally may be quite difficult 
regarding the gluons. For the sea quarks, however, there may be a chance to have further 
constraints from the NA50 data for Drell-Yan in Pb-Pb collisions at = 17.2 GeV at the 
CERN-SPS, provided that the precision in the invariant mass distribution will be sufficient 
at large masses. If the sea quarks show an EMC effect, this should serve as an additional 
constraint for the gluons through the evolution. 

The region at x < 0.01 involved assumptions of saturation of shadowing, and at very small 
values of x the gluons were expected to have the same shadowing as (or the sea quarks). 
To pin down the nuclear gluon shadowing at small values of x, high precision data on the 

dependence of F2 / F2 , F2 I ^2 would be needed especially at 0.01, where the gluons 
start to determine the evolution of the sea quarks more strongly |T^, In the future. 



perhaps also the gluon initiated processes like open charm and J/\l/ production in DIS |^ 



and in pA collisions |3^, could help in constraining the gluon antishadowing peak further. 

Related to the small-x physics, the role of higher twist effects like the GLRMQ fusion 
corrections needs still to be studied further together with the results from HERA. Also 
the NLO analysis, both in the cross sections and in the evolution, should be performed. Most 
importantly, however, our analysis would benefit from an inclusion of a more quantitative 
estimates of uncertainties and propagation of errors. 
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Figure 1: The nuclear ratios Rf{x, Q^) for individual parton flavours i = g, uv,u, d, s, c of a lead 
nucleus A = 208 as functions of x at fixed values of Q"^ = Ql = 2.25 GeV^ and = 10000 GeV^ as 
obtained by using the GRV-LO |25] distributions (solid lines) and the CTEQ4L [^6| distributions 
(dotted lines) for the free proton. The dashed lines show our numerical parametrization (EKS) 
to the nuclear effects obtained in the GRV-LO case. The ratios Rf^ are almost identical to R^y: 
and are not shown. The charm ratios are presented only for = 10000 GeV^, since the charm 
distributions are generated only above our Ql The ratios R{^ at = 10000 GeV^ behave as R^, 
and are not shown. For practical purposes the set-dependence of the ratios Rf{x, Q^) is negligible. 
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Figure 2: The nuclear ratios Rf{x,Q^) for individual parton flavours i = g,uv,u,d,s,c as 
functions of the mass number A at a fixed value of = 10 GeV^ and at fixed values of x = 
10"'^, 10~^, 10~^, 10~^ (see the panel for the strange quarks). The notation of the curves is the 
same as in Fig. || and the markers show the nuclei for which we have made the computation. The 
set-dependence is the largest for large nuclei but keeps within ~5 %. 



12 



c 



1.1 
1.05 

1.0 
0.95 

0.9 
0.85 

0.8 

1.1 
1.05 

1.0 
0.95 

0.9 
0.85 

0.8 

1.1 
1.05 

1.0 
0.95 

0.9 
0.85 

0.8 

1.1 
1.05 

1.0 
0.95 

0.9 
0.85 

0.8 

1.1 
1.05 

1.0 
0.95 

0.9 
0.85 

0.8 
0.75 



~i — I I 1 1 1 1 II — I — I I 1 1 1 III 




1 



x=0.0125 



I I I 1 1 



I I I I N il 




x=0.035 



I I I 1 1 



I I I I M il 



.x=0.070 




x=0.175 



I I I 1 1 



I I I I N il 



x=0.45 

I 



iL 



I I 1 1 1 1 1 — I — I I 1 1 1 III — I — I — I 1 1 1 1 II — I — I I 1 1 1 1 1 



x=0.0175 



I I 1 1 n il 




x=0.045 



I I I I III 



.T=0.090 



I I I I III 



x=0.25 



I I I UN 



x=0.55 



x=0.025 




.T=0.055 




x=0.125 



x=0.35 



x=0.10 

I 



10 



100 1 



10 



100 1 



10 



100 



1.1 

1.05 
1.0 
0.95 
0.9 
0.85 
0.8 
1.1 
1.05 
1.0 
0.95 
0.9 
0.85 
0.8 
1.1 
1.05 
1.0 
0.95 
0.9 
0.85 
0.8 
1.1 
1.05 
1.0 
H 0.95 
0.9 
H 0.85 
0.8 
1.1 
1.05 
1.0 
0.95 
0.9 
0.85 
0.8 
0.75 



Figure 3: The ratio / F2 as a function of at several different fixed values of x. The 
data is from Q. The figure contains four calculated curves: two of the curves correspond to the 
"exact" results obtained with the nuclear ratios for the GRV-LO distributions (as in [l^) and for 
the CTEQ4L-distributions separately, and two curves correspond to the results obtained with the 
GRV-LO and CTEQ4L distributions multiplied by our numerical parametrization of Rf{x,Q'^). 
There is no significant difference between the calculated curves. 
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Figure 4: The ratio of differential Drell-Yan cross sections in pA and pD from Eq. as a function 
of X = X2 for i|C, IgCa, |^Fe and if|W. The open squares show the E772-data [||], and the filled 
symbols stand for the calculated ratios i?^y(x,(5^) at {x,Q^) corresponding to the experimental 
values ||3^. The circles show R^y computed with the nuclear ratios obtained separately for the 
GRV-LO set (big circles) and for the CTEQ4L set (small circles). The results obtained by using our 
numerical parametrization (EKS) of Rf together with the sets GRV-LO and CTEQ4L are shown 
by triangles and diamonds, correspondingly. As seen from the panel for tungsten, the differencies 
between the two parton distribution sets used for the free proton are larger than the error from 
using the set-independent parametrization for the nuclear effects Rf. 
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